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1.  Introduction 


Complex  ABCD  ray  matrix  techniques  combined  with  the  paraxial  approximation 
of  the  Huygens-Fretnei  formulation  of  wave  optics  is  a  useful  tool  In 
analyzing  analytically  the  performance  of  a  rather  general  dees  of  optical 
systems^^^  This  method  expresses  the  complex  optical  field  in  the 
observation  (or  output)  plane  as  an  Integral  over  an  initial  (or  Input) 
plane  of  the  given  source  distribution  multiplied  by  a  propagation  kernel, 
which  iu  a  function  of  the  ABCD  matrix  elements  of  the  complete  optical 
system  between  the  input  and  output  plane  .and  propi^ties  of  the  Intervening 
medium  which  the  light  propagates  through.  Thus,  for  example,  the  effects  of 
a  random  medium  (e.g.,  turbulence)  and  incoherent  sources  can  be  described 

(2  3) 

by  using  this  technique  .By  employing  such  methods  one  is  able  in  a 
straight  forward  manner  to  obtain  expressions  for  all  the  statistical 

correlation  functions  that  pertain  to  the  optical  system,  source 
distribution  and  intervening  media  under  consideration. 

The  use  of  complex  ABCD  matrix  elements  permits  modelling  of  limiting 
apertures  (e.g.,  thin  lens,  field  stops,  and  finite  sized  measurement 
apertures)  as  soft  (i.e.,  Gaussian)  apertures.  In  many  applications  (e.g., 

imaging  and  Fourier  transform  systems)  such  modelling  yields  useful 

engineering  analytical  approximations  to  system  performance  that  can  be  used 
for  designing,  sizing  and  sealing  of  optical  systems.  In  addition,  the 
results  can  generally  be  cast  in  a  physically  intuitive  form,  where,  for 
example,  the  effects  of  geometric  optics  and  diffraction  are  readily 

apparent. 
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In  this  psper  we  consider  homogeneous  media  and  derive  in  sec.  2  expressions 
for  the  mean,  variance  and  corresponding  correiation  function  of  the 
irradiance  distribution  resulting  from  an  incoherent  source  that  has 
propagated  through  a  complex  (axially  symmetric)  ABCD  optical  system.  In 
particular,  we  obtain  a  general  expression  for  the  mutual  coherence  function 
which  can  be  regarded  as  a  generalization  of  the  van  Cittert-Zemlke  theorem 
to  complex  paraxial  ABCD  systems.  We  also  derive  an  expression  that  relates 
the  number  of  speckle  correlation  cells  contained  within  a  measurement  area 
to  the  parameters  of  the  ABCD  system.  In  partlcuiar,  an  expression  for  the 
maximum  number  of  independent  intensity  measurements  that  the  optical  system 
allows  is  obtained. 

In  sec.  3  we  consider  the  effects  of  wave  front  decorrelation  on  fringe 
formation  in  speckle  interferometric  systems.  Such  interferometers  brings 
about  a  real-time  technique  for  non-destructive  measurements  of  deformations 
and  contours  of  objects  that  scatter  laser  light  diffusslvely.  To  make 
measurements  using  double  exposure  speckle  interferometry,  a  speckle-pattern 
exposure  of  the  object  is  taken  in  one  position.  The  object  is  then  deformed 
(i.e.,  put  under  load)  and  another  exposure  is  taken.  These  two  speckle 
patterns  are  subtracted,  and  their  difference  is  squared  to  obtain  speckle 
correlation  fringes  that  correspond  to  the  objects  movement  and 

(4) 

deformation  .  Here  we  examine  the  effects  of  load  induced  object  tilt 
around  an  In-plane  axis,  in-plane  and  out  of  plane  object  displacement  on 
decorrelation  of  the  Interferograms  operating  in  the  subtraction  mode. 

Decorrelation  and  the  corresponding  fringe  visibility  in  image-plauie 

electronic  speckle-pattern  correiation  interferometers  has  been  considered 

(5) 

recently  by  Owner-Petersen  .  However,  it  can  be  shown  (see  Appendix)  that 
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the  correlation  function  derived  in  Ref.-  5  is  valid  only  for  the  cate  where 
the  number  of  independent  speckle  modes  N  is  much  larger  that  unity  (i.e.,  N 
•*  tt).  Here  we  derive  expression  for  the  correlation  fxmetion  for  both  an 
arbitrary  ABCD  system  and  finite  N.  In  particular  it  is  shown  that  as  N  *»  1, 
the  correlation  coefficient  is  equal  to  unity.  That  is,  in  this  limit  the 
performance  of  phase  difference  speckle  interferometers  do  not  degrade  under 
load  induced  object  deformations. 

2.  Intensity  Statistics  for  Incoherent  Sources 

We  consider  a  completely  incoherent,  finite  sized,  quasi-monochromatic, 
polarized  optical  source  located  in  an  input  plane  and  assume  that  a 
paraxial  optical  system,  characterized  by  a  complex  ABCD  ray  matrix, 
occupies  the  space  between  the  primary  source  and  an  observation  plane.  For 
definiteness,  we  consider  an  axially  symmetric  paraxial  optical  system  that 
Includes  a  radially  Increasing  optical  amplitude  loss  with  respect  to  the 
optic  axis.  In  particular,  we  consider  a  general  paraxial  system  which 
Includes  a  Gaussian  variation  in  amplitude  transmission  across  the  axis.  It 
has  been  shown  that  such  systems  can  be  described  by  a  complex  valued  ABCD 
ray  matrix 

By  using  a  Huygens-Fresnel  integral  formulation,  whose  propagation  kernel  is 
a  function  of  the  system  ABCD  matrix  elements,  one  can  express  the  complex 
optical  field  in  the  observation  plane  of  an  axially  symmetric  system  where 
the  input  and  output  planes  are  in  the  same  medium,  as 

r  T  r  * 

U(s)  -  I-  ejcpf-lJtL;Jd>  U^(r) 
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X  exp  -  2r-p  ♦  Ar^ 


(1) 


where  k  it  the  opticel  wave  number  in  free  apace,  L  ia  the  optical  diatance 

along  the  s^axia,  r  and  g  are  coordinate  veetora  tranaverae  to  the  z-axia  in 

the  aource  and  obaervatlon  plane,  reapeetively,  t  ia  the  on-axia 

tranimiaaion  amplitude  factor  and  A,B,  and  D  are  the  overall  complex  ray 

matrix  elementa  for  the  ayatem  under  conaideration  from  the  input  to  the 

output  plane.  Without  loia  of  generality  we  can  aaaume  that  t  ■  2. 

0 

We  now  apply  Cq.  (1)  to  the  epeclal  caae  where  the  aource  ia  incoherent  and 
obtain  general  expreaaiona  for  the  mean  and  second  moment  of  the  obaerved 
intenaity.  In  any  experimental  meaaurement  of  the  intenalty  one  employe 
detectora  of  finite  aize.  Following  Goodman  the  inatantaneoua  meaaured 
intensity  is  expresaed  aa 


•  5  Jd*p  W(s)  1(2).  (2) 

where  p  ■  (p  ,  p  ),  Jfg;  ■  It/fg.)!*,  W(_p)  ia  a  real  and  poaltlve  meaaurement 
•*  ¥ 

aperture  weighting  function,  and  S  ■  Ju^CgM^p.  Due  to  the  preaence  of  W(2) 
the  integration  Indicated  in  the  above  can  be  eonaidered  to  be  over  the 
entire  observation  plane. 

From  Eqa.  (1)  and  (2)  the  mean  meaaured  intenaity  can  be  expressed  as 

<I^>  •  I  Jd*p  W(2)  <I(2»  (3) 
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where 


Jd*r  Ij!r)  §xp  [~2k  Im(l/B)£‘rJ 

X  exp[k  Im(A/B)r^],  (4) 

where  Im  denotes  the  Imagintry  pert,  and  the  angular  brackets  Indicate  on 
ensemble  average  over  the  statistics  of  the  source.  In  obtaining  Eq.  (4)  we 
used  the  relation,  valid  for  an  incoherent  source  radiating  into  a 
half-space,  (k^/2n><U  (r}U*(rJ>  -  I  [(r*r)/2J8(r-rJ,  where  8(r)  is  the 
two-dimensional  Dirac  delta  function  and  ir  the  radiated  power  per  unit 
area  of  the  source.  By  considering  the  field  of  a  point  source  located  at  a 
point  r  in  the  input  plane.  It  can  be  shown  for  a  system  with  no  gain  and  a 
radially  increasing  loss  across  the  optic  axis  that  both  Im(A/B)  and  Im(D/B) 
are  lessn4Xthan  or  equal  to  zero.  As  a  result,  the  Integrals  In  Eqs.  (3)  and 
(4)  are  finite  for  any  practical  source  distribution  of  Interest. 

In  a  similar'  manner,  the  mean  square  of  the  measured  intensity  can  be 
expressed  as 

•  -2  fS) 

Assuming  circular  complex  Gaussian  statistics  for  the  underlying  source 
field,  it  can  be  shown  that  the  intensity  autocorrelation  function  can  be 

.xpr....d  » 
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where 


(7) 


if  the  abiolute  magnitude  of  the  complex  degree  of  coherence.  The  quantity  T 
it  the  mutual  coherence  function  of  the  field  in  the  obaerration  plane  and  ie 
given  by 


-  <U(s^)  U  (s^)> 


1  1* 

(Jc 

(D  2 

f 

■ 

1  255  1  “f 

■  1 

■ 

UPi* 

?”») 

■ 

(r) 


X  exp 


Ikr* 


e.  J.i 

r-pl 


(8) 


Equation  8  le  the  generalization  of  the  mutual  coherence  function  for 
propagation  in  a  homogeneoua  medium,  given  by  Eq.  (20)  of  Ref.  7  chapter 
10.4.2  ,  to  complex  ABCD  optical  eyatema.  Thua,  Eq.  (8)  can  be  regarded  aa  a 
generalization  of  the  van  Cittert-Zernike  theorem  to  complex  paraxial  ABCD 
optical  aystems. 


Examination  of  Eq.  (8)  reveala  that  the  atatiatlca  of  intenaity  for  a  complex 
ABCD  ayatem  are  not  atatienary  but  rather  depend  on  abaolute  poaltion  and 
gj*  For  a  given  circularly  aymmetric  ABCD  ayatem  containing  aoft  aperturea 
the  mean  irradiance  and  mutual  coherence  function  la  given  by  Eqe.  (4)  and 
(8),  reapectively.  The  correapondlng  meaaured  mean  and  mean  aquare  Intenaity 
are  given  by  Eqa.  (3)  and  (5).  Ehie  to  the  aaaumptlon  of  an  underlying 
Cauialan  atatiatlca  for  the  complex  field  all  higher  momenta  of  the  intenaity 
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can  be  expressed .  by  suitable  products  of  the  corresponding  first  and  second 


moments. 


Another  quantity  of  interest  is  the  number  of  speckle  correlation  cells  N 
within  the  measurement  aperture.  Following  Goodman^^^  this  quantity  Is  given 
by 


"■j'ZTy 


> 

M 


<I*>  “  <l  >* 

M  M 


(9) 


Substituting  Eqs.  (3)  and  (5)  into  Cq.  (9)  yields 

tSW(z)  <I(2»d*p]* 

N  >  . - . . . . - . . . . 

ssw(2^)w(z^)<i  ( j  ;><K  I  |  ^d\d\ 


(10) 


where  </('g.)>  and  era  given  by  Eqs.  (4)  and  (7),  respectively. 

In  order  to  proceed  we  specialize  to  the  ease  where  the  intensity 
distribution  of  the  source  is  of  gausslan  shape  given  by  l(r)  ■ 
I  txp(-2r^/r^^).  In  addition  we  consider  the  situation  where  the  beam  spot 
size  in  the  observation  plane  Is  much  smaller  than  the  extent  of  the 
measurement  aperture;  that  is.  we  assume  W  ■  constant  for  K^)  *  0,  Note,  if 
the  measurement  aperture  is  characterized  by  a  Gaussian  shaped  weighting 
function,  it  can  be  formally  included  directly  In  the  ABCD  matrix  elements  as 
being  the  last  optical  element,  immediately  In  front  of  the  observation 
plane. 
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Ui?on  performing  the  Integration  indicated  in  Eq.  (8)  it  can  be  shown  that  the 
mutual  coherence  function  is  given  by 


r  'vv*!  r  f'!*  1 

rcp.p;-  — ^  - ! i.  exp  — ^ 

*  *  JtV  ‘  p*  «* 


(11) 


exp 


[‘  - 1 


where 


s* .  iiii! .  j 

0  Jc*p*  ^ 


(12) 


and 


c*  -p* 
0 


Ifii 


k  r* 


(13) 


The  corresponding  expression  for  the  beam  intensity  profile  and  the  magnitude 
of  the  degree  of  coherence  are  given  by 


■  -j-l  exp  2p*/M*j 


(14) 


and 


|r|  -  expj-  <2-  g//p*] 


(IS) 
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The  quantity  p  is  a  measure  of  the  lateral  extent  of  the  region  in  the 
0 

observation  plane  where  the  optical  field  it  correlated  (and  is  hence  a 
measure  of  the  mean  speckle  size)  while  the  quantity  u  is  a  measure  of  the 
beam  intensity  spot  size.  For  a  given  ABCD  system,  these  quantities  can  be 
obtained  directly  from  Eqs.  (12)  and  (13). 

The  number  of  speckle  correlation  cells  within  the  beam,  as  obtained  from  Eq. 
(10),  is  given  by 


N  m  z!L 

ol 


(16) 


As  an  illustrative  ABCD  system  we  consider  an  optical  system  that  images  a 
Gaussian  shaped  incoherent  source  onto  an  observation  plane.  A  positive  thin 
lens  of  Gaussian  transmission  radius  o  is  located  a  distance  to  the  right 
of  the  source.  At  a  further  distance  a  Gaussian  shaped  aperture  (l.e.,  a 
measurement  aperture)  of  radius  o-^  is  Immediately  to  the  left  of  the  image 
plane.  The  appropriate  ray  matrix  elements  can  be  calculated  in  a 
stralghforward  manner  with  the  result  that 


A  D 


(17a) 


B  • 


2t 


and 


D  - 


2t  L 


k<r‘ 


(17b) 


(17c) 
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Substituting  £qs.  (17)  into  Eqs.  (12)  an'l  (13)  yields  that 


and 


where 


tiL 

kV 


41 


1  * 


t 


Jc*rV 

■ 


r/fj 

2 

2 

P*  ♦ 

■ 

2 

TtoT 

k  1 

1 

(18) 


(18) 


(20) 


is  the  spot  size  in  the  observation  plane  (i.e.,  image  plane)  for  an  Infinite 
measurement  aperture  radius.  Examination  of  Eq.  (20)  reveals  that  there  are 
two  contributions  to  the  first  term  under  is  a  geometric  magnification  of 
the  source  size  and  the  second  term  represents  the  effects  of  diffraction 
from  the  finite  size  of  the  imaging  lens. 


Equation  (19)  gives  the  effect  of  a  finite  measurement  aperture  on  the  image 
size.  As  expected,  for  values  of  much  larger  (smaller)  than  the 
resulting  image  size  is  given  by  u  (e* ).  On  the  other  hand  the  coherence 

V  A 

length  is  independent  of  e*^,  as  it  should. 

The  number  of  correlation  cells  contained  in  the  imaged  spot  is  obtained  from 
Eqs.  (16],  (18),  and  (19)  as 
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AT  ■  1  + 


(21) 


2  (v  /p  f' 

4  O 

1  ♦  /w 

•  CD 

and  it  plotted  in  Fig.  1  at  a  function  of  (!•«••  the  ratio  of  the 

meaturement  aperture  to  the  coherence  length)  for  varloue  faluet  of  e  /u  . 

For  fixed  the  number  of  correlation  oella  inereaae  for  deereaalng  ealuee 

of  p  ,  the  corresponding  increase  being  .progressively  less  for  increasing 
0 

values  of  o*  /«  . 

4  ip 

As  a  further  illustrative  example,  we  consider  a  Fourier  transform  geometry 
with  a  finite  sized  thin  transform  lens  of  ‘‘positive  focal  length  f  and 
gausslan  radius  c.  As  in  the  previous  example  a  gausslan  shaped  measurement 
aperture  of  radius  it  assumed  to  be  Just  in  front  of  the  Fourier  plane.  In 
this  ease  the  relevant  matrix  elements  are  given  by 

(22a) 

(22b) 

(22c) 

and  u^,  the  spot  size  in  the  Fourier  plane,  is  given  by  Eq.  (19)  where 


..ill 


k9-* 


.f-in 

hr* 


B  •  f 


and 


4r  .  £0- 


acre  r 


L^L 

e*  e» 


Substituting  Eqt.  (22)  into  Eqt.  (12)  and  (13)  yields  that 
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Figure  1.  Hw  number  of  Intensity  correlation  cells  for  Imaging  and 
Fourier  transform  geometries  contained  in  the  measurement  aperture 
as  Ainction  of  the  ‘^aperture  radius:”  (normalized  to  the  lateral  coher¬ 
ence  length  In  the  observation  plane)  for  various  values  of  Oa/co*.. 
For  the  Imaging  and  Fourier  transform  geometries,  the  quantities  Po 
and  (Dm  are  given  by  Eqs.  (18)  and  (19),  and  Eqs.  (23)  and  (24), 
respectively. 
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(24) 


The  number  of  correlation  cells  contained  in  the  measurement  aperture  is 
formally  the  same  as  Eq.  (21),  with  p  and  u  given  by  Eqs.  (23)  and  (24), 

O  00 

respectively. 


Finally,  consider  the  case  of  an  extended  incoherent  source  (i.e,,  r  ->  m).  Of 

O 

particular  interest  is  the  number  of  correlation  cells  N.  For  r  -»  «,  N  is 

O 

independent  of  the  source  and  depends  on  the  wavelength  and  properties  of  the 
optical  system  only.  Since  the  underlying  probability  density  for  the  complex 
field  has  been  assumed  Gaussian,  N  also  gives  the  number  of  statistically 
independent  speckle  modes.  For  ->  ob,  N  can  be  interpreted  as  the  n^imum 
number  of  independent  intensity  measurements  that  the  optical  system  allows. 
For  an  arbitrary  complex  ABCD  system  the  maximum  number  of  independent 
intensity  measurements  allowable  is  given  by 

N  =  1  +  - - — - -  (25) 

Jm(BA  )Im(BD  )  -  (ImB) 


Note,  for  a  point  source  (i.e,,  r^-)  0)  the  corresponding  number  of  modes,  as 
expected,  is  equal  to  unity.  For  the  imaging  smd  Fourier  transform  geometries 
N  is  given  by 

mAx 


N 

max 


fko*  (T 


I  + 


2d 


(25) 


where  d  =  and  f  for  the  imaging  and  Fourier  transform  geometry, 
respectively. 
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3.  Correlation,  Functions  in  Phase  Shifting  Interferometers 

A  basic  characteristic  of  phase  shifting  interferometer  operating  in  the 

subtraction  mode,  being  electronic  or  holographic,  is  the  formation  of 

secondary  fringes  by  taking  the  difference  and  squaring  the  two  primary 

(4) 

interferograms,  before  and  after  loading  .  Then  the  ensemble  averaged 
secondary  fringe  pattern  is  obtained  as 

<(I  -  I  )^>  =  <lS  +  <lS  -  2<I  I  >,  (26) 

»  b  *  b  kb 

where,  for  notational  convenience,  we  have  suppressed  the  functional 
dependence  on  observation  plane  coordinates  p,  the  subscripts  b  and  a  denote 
before  and  after  loading,  and  angular  brackets  denote  the  ensemblef  average 
(i.e.,  over  all  possible  realizations  of  the  source  distribution  of  scattered 

light).  Now  each  of  the  interferograms  I  and  I  ,  found  from  two  primary  wave 

A  b 

fronts  and  U^,  is  given  by 

I  =  |U+  U  1^  =  1+  I  +  2Re[UU  *]  (27) 

'12'  12  12 

where  I  =  |U  |^.  Assuming  that  U  and  U  are  stochastically  independent 

1,2  1,2  1  2 

and  that  the  intensity  of  a  speckle  pattern  obeys  negative  exponential 

statistics,  M.  Owner-Petersen  has  calculated  the  secondary  fringe  pattern  and 
corresponding  visibility  V  of  the  correlation  fringes  for  the  cases  where 
is  speckled  and  may  be  speckled.  The  situation  where  and  are 
speckled  and  decorrelate  under  load  is  relevant  both  to  the  in-plane 

displacement  sensitive  interferometer  and  to  the  out-of -plane  gradient 
interferometer,  where  both  the  primary  wavefronts  and  are  scattered 
from  the  surface  under  test.  The  situation  when  only  decorrelates  under 
load  is  relevant  to  the  c.’t-of-plane  displacement  sensitive  interferometer, 
where  represents  a  deterministic  reference  wave.  Central  to  the  quality  of 
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Preceding  Page  Blank 


the  eecondary  fringei  in  \':he  verioue  Interferometen  is  the  correlstion 
coefficient  of  the  complex  fields  in  the  observation  plane  before  and  after 
loading.  This  quantity  is  defined  as 

<UA}*> 

^  / - T-  / - r  <28) 

where  it  is  assumed  that  both  U.  and  U.  have  the  same  eorrelation 
funetlon<^\  In  general  t  ie  eompiex  with  0«|3r|ai,  Thus  the  effects  of 
decorrelation  are  a  distortion  and  corresponding  reduction  of  fringe 
visibility  in  the  observation  plane  (e.g.,  V  o  as  |r|4  o).  Owner-Petersen 
derives  an  explicit  expression  for  the  decorrelation  coefficient  in  a  direct 
image-collecting  optical  system.  However,  as  indicated  in  the  appendix,  the 
results  given  in  Ref.  5  apply,  strictly  speaking,  to.  the  case  where  the 
number  of  modes  collected  by  the  imaging  system  are  very  large. 

Specifically,  r  is  the  cross-correlation  function  of  the  complex  optical 
field  before  and  after  loading  observed  at  a  point  in  the  observation  plane, 
l.e.,  in  practice  at  one  resolution  or  pixel  element  in  a  detector  array. 
Here  we  compute  r  for  three  types  of  deformation,  which  for  sufficiently 
small  loadings  will  be  Independent:  (1)  tilt  about  an  in-plane  axis,  (2) 
in-plane  translations,  and  (3)  a  translation  parallel  to  the  optic  axis.  Here 
we  assume  an  object  is  illuminated  by  a  laser  beam  of  mean  waveveotor  ]^.  The 
scattered  light  of  wavevector  h  (dong  the  optic  axis  propagates  through  a 
complex  ABCD  and  the  irradianee  is  recorded  by  a  suitable  detector  array  in 
the  observation  plane.  It  is  assumed  that  the  object  is  sufficiently  rough 
such  that  the  reflected  light  can  be  assumed  to  be  spatialiy  completely 
incoherent  in  the  object  plans.  Specifically,  we  consider  that  the  reflected 
light  is  given  by  u(r)  •  p(r)  txp  [i(k  -  k  )•  zJ,  where  p(r)  is  the  (complex) 
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reflection  coefficient  of  the  object.  For  en  incoherent  object  radiating  Into 
a  half  space  we  have  that  <p(r  )p*(r.)>  ■  (2ii/k*)  d(r  -r  )I  ((r  +r  )/2j, 
where  k  is  the  magnitude  of  the  wavenumber,  1^  is  proportional  to  the 
brightness  distribution  of  the  source  and  ia  assumed  here  to  be  of  the  form 
expl-2rVr*l,  where  r^  is  the  source  spot  radius. 

3.1  In-Pling  .Tilt 

In  this  situation,  the  reflected  optical  field  after  loading  in  an  object 
plane  transverse  to  the  r-sxis  is  given  by  (assuming  small  tilt  angles  9) 

U  (r)  -  expIIk^erlUJ^),  (29 

where  U.  (c)  la  the  reflected  field  before  loading, 

9  ■  Cl  +  cosa),  (30) 

and  a  is  the  angle  between  Jc^  and  the  9-axis  (i.e.,  the  angle  of  incidence  of 
the  illumination  beam). 

The  correlation  function  is  obtained  by  employing  Eq.  (29)  in  Eq.  (1)  and 
substituting  the  resulting  expressions  into  Eq.  (28).  The  ensemble  average  is 
readilly  obtained  by  noting  that  <U.(r  )U/(r,)>  ■  (2«/k*)I  I(r  ♦r,)/2] 

Assuming  a  Gaussian  shaped  source  spot  size,  the  resulting  Gaussian  integral 
can  be  integrated  analyticaliy.  After  a  straight  forward  calcuiation,  it  can 
be  shown  that  (assuming  a  tilt  about  the  y-axis) 

r  -  Irje  (31) 
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where 

ly|  ■  expl-  eVe*i, 


(32) 


6 


2|B0 


(33) 


« 


t 


4BIm(B)p^9 


(34) 


is  given  by  Eq.  (12),  B  is  the  B-mstrix  element  of  the  optical  system 

O 

before  loading,  and  p^  is  the  projection  of  the  observation  coordinate  vector 

onto  the  x-axis.  Equations  (31)-(34)  give  the  load  induced  object  tilt 

correlation  coefficient  for  a  complex  ABCD  optical  system  applicable  to 

phase-shifting  speclcle  interferometry.  Note,  for  a  point  source  (i.e.,  r^<4  0, 

p  4  «)  we  obtain,  as  expected  intuitively,  that  y  4  1. 

0 


3.2  In-Plane  Translation 


In  this  situation  the  load  induces  a  uniform  displacement  of  the  object  by  an 
amount  Ar  perpendicular  to  the  B-axis.  Thus,  the  reflected  object  field 
after  loading  can  be  expressed  as 

U^(r)  -  U^(r-Ar)  (35) 

By  following  the  same  procedure  used  in  see  3.1,  it  can  be  shown  that  the 
in-plane  translation  induced  correlation  coefficient  is  given  by 

y  ■  |y|e‘*expll(k-k  )  •  Ar2.  (36) 
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where 


»  ■  exp 


(Ar)'* 

.a 


(37) 


•  ■  At  +  Cj^  E  *  ^ 


Ar 


4Re(BA») 

kr*p* 


(38) 


(39) 


(40) 


and 


c  -  )c  Re(l/B)  - 
a 


4Im(l/B)Re(BA*) 


(41) 


Equation  (35)-(41)  gives  the  correlation  coefficient  for  load  Induced 
ln*plane  translation  for  complex  ABCD  eystems.  For  a  point  eouree  (r^-^  0)  on 
the  optic  axis  we  obtain  that  |r|  ■  1.  ■  k  Re(A/B)/2,  and  ■  h  Re(I/B), 

which  is  identical  to  what  one  would  obtain  from  Eq.  (1)  directly.  As 
expected  physically,  the  correlation  coefficient  for  a  point  source  is  a  pure 
phasor. 


ia.,TriMllttgni  BtTilltl  tB  \tit  Mfll 


When  the  load  induces  a  translation  along  the  z-axis  by  an  amount  Az,  the 
field  after  loading  can  be  obtained  from  the  oorresponding  field  before 
loading  by  replacing  the  ABCD  matrix  elements  by  A,  B,  S,  and  6,  where 
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A  B  1  _  [a  b1[i  Azl 

c  6  J  ‘  [cdJLo  1  J 

That  It.  it  given  by  Eq.  (1)  with  the  matrix  elementt  (A  i  £  6).  where  X  ■ 
A,  I  ■  B^AAx,  C  ■  C,  and  5  «  D  ^  CAz. 


Simllarliy,  by  following  the  tame  procedure  outlined  In  eeo  3.1.  It  can  be 
thown  that  the  correlation  coefficient  for  tranalatione  along  the  optic  axil 
It  given  by 


y  -  a  exp  (Kk-k  )  Azl  expllklRe(D/B)  -  Re(D/S)lp*/2)  expt-c  pWl.  (42) 

X  w 


where  u,  the  tpot  tlze  In  the  obtervation  plane  before  loading,  it  given  by 

Eq.  (13). 


a  ■ 


i^!i  ’ 


(43) 


p: 


-  4  2Az 


4AB* 


(krj* 


1^ 

k 


(44) 


p  it  given  by  Eq.  12,  p  ^  it  obtained  from  Eq.  (12)  by  replacing  A  and  6  by 
A  and  §,  retpectlvely,  and  it  a  (complex)  conttant  with  |e^|  <1,  and 
Re(Cj)  a  0.  Although  a  rather  lengthy  expreitlon  for  at  a  function  of  Az 
and  A...  D  can  be  written  down  it  it  of  no  great  concern  here.  Thit  can  be 
teen  by  recognizing  that,  in  practice,  meaningful  information  from  the 
eecondary  fringe  are  obtained  In  the  regime  of  the  obtervation  plane  where 
the  recorded  inteneity  it  relatively  high  (i.e..  for  p  <  w). .  Rettrlctlng  our 
attention  to  valuet  of  p  near  the  center  of  the  tpot  (i.e.,  near  p  •  0),  the 
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magnitude  of  the  correlation  coefficient  for  displacements  parallel  to  the 
optic  axis  is  equal  to  the  amplitude  coefficient  a,  given  by  Eq.  (43).  For  a 
polntsource.  a  4  1,  the  correlation  coefficient  is  a  pure  phasor. 

3.4  AppUcitioni 

We  now  apply  the  results  of  secs  3. 1*3.3  to  three  optical  systems  of  concern 
to  phase  difference  speckle  interferometry  and  obtain  analytic  expression  for 
the  magnitude  of  the  correlation  coefficient  |  r  |  •  This  quantity  relates 
directly  to  fringe  visibility  see  Eqs.  (15)  and  (17)  of  Ref.  5). 

As  Indicated  schematically  in  Fig.  2,  the  clean  imaging  optical  system 
consists  nf  a  large  positive  thin  lens  of  focal  length  f^  placed  a  distance 
f^  to  the  right  of  a  Gaussian  shaped  source  of  radius  r^.  At  a  distance  f^  to 
the  right  of  this  lens  is  a  Gaussian  shaped  limiting  aperture  of  radius  o'.  At 
a  further  distance  f^  to  the  right  of  this  aperture  is  a  large  positive  thin 
lens  of  focal  length  f^.  Finally,  the  observation  plane  is  located  a  further 
distance  f^  to  the  right  of  this  lens.  An  imaging  system  of  this  type  is 
convenient  since  by  varying  the  limiting  aperture  radius  o*.  one  can  match  the 
target  resolution  to  ensure  a  fully  resolved  speckle  pattern  across  the 

detector  array.  Furthermore,  since  the  limiting  aperture  is  placed  in  the 

Fourier  plane  the  clean  imaging  system  is  relatively  insensitive  to 

vignetting  effects,  and  quadratic  phase  factors,  present  in  direct  imaging 
systems,  are  absent  here  leading  to  mathematical  simplifications  of  the 

complex  field  in  the  observation  plane. 
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The  appropriate  ray  matrix  elements  are  given  by;  A  ■  f  /f  ,  B  ■  21f  f  /Iw*, 

2  1  12 

and  D  «  ^/^2'  inean  speckle  size  and  magnitude  of  the  spot  size  u  are 
given  by 


where  N,  the  number  of  modes  "captured"  by  the  system,  takes  the  form 


N  -  1  + 


-Wf 


(47) 


Next  consider  a  direct  imaging  system  where  a  positive  lens  of  focal  length  f 
and  Gaussian  radius  c*  is  located  a  distance  to  the  right  of  a  source  of 
radius  r  .  At  a  further  distance  L,,  to  the  right  of  the  lens  (where  L’^ 
>  f*  )  is  the  observation  plane.  The  appropriate  ABCD  matrix  elements  are 

given  by  Eqs.  (17)  for  o*  ®.  It  is  easily  seen  that  the  quantities  p  ,  w, 

®  0 

and  N  are  identical  to  those  obtained  above  for  the  "clean"  imaging  system 

Finally,  v;e  consider  a  Fourier  transform  optical  system.  Such  a  system  may  be 
of  interest  to  measurement  of  surface  rotation  ana  tllt^^^  Tiziani  shows 
that  under  these  conditions  the  rigid  body  tilt  is  measured  independently  of 
any  lateral  displacement.  The  appropriate  ray  matrix  elements  are  given  by 

I 

Eqs.  (22)  and  we  readily  obtain  from  Eqs.  (23)  and  (24)  that 
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and 


(49) 


where  the  number  of  modes  N  it  given  by  Eq.  (47). 

In  Table  1  we  present  the  corresponding  results  for  the  correlation 
coefficients.  For  single  mode  systems  (i.e.>  N  •*  1,  p  ■»  •)  the  correlation 
coefficient  is  unity  and  fringe  visibility  will  be  unaffected  by  lead  Induced 
deformations.  On  the  other  hand,  for  direct  imaging,  and  N  »  1,  the  results 
presented  in  Table  1  for  both  the  in-plane  tilt  and  displacement  agree 
qualitatively  with  the  results  of  Ref.  5  for  the  case  where  the  imaging  lens 
is  in  the  limit  of  Fresnel  diffraction  (l.e.,  «  kff*). 

To  illustrate  the  effects  on  decorrelation  due  to  finite  N  we  consider  clean 
imaging  and  plot  the  |r|  as  a  function  of  both  9  and  Ar,  normalized  to  the 
respective  decorrelation  parameter  for  N  ■»  «.  Since  both  9  and  r .  depend  on 

1  d 

N  in  an  Identical  manner  we  obtain  a  single  family  of  curves,  which  are 
plotted  in  Fig.  3.  Examination  of  Fig.  3  indicates  that  for  Naio,  one  may 
employ  the  decorrelation  lengths  for  both  tilt  and  in-plane  displacement 
obtained  for  N  -»  «  with  negligible  error. 

For  local  induced  deformations  along  the  optic  axis,  we  see  from  Table  1  that 
|jr|  decorrelates  with  the  characteristic  length  z^  (where  z^  is  of  the  order 
the  depth  of  field)  for  the  clean  and  direct  imaging  (for  L^«  Im^)  and  with 
the  fresnel  length  for  the  Fourier  transform  system. 
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Figure  3.  The  magnitude  of  the  conelatlon  coefficient  for  Imaging 
systems  as  Amotion  of  normalized  deformations  for  various  values 
of  N.  For  In-plane  tilt  and  ln<plane  displacement  x  ■  9/6tao  and 
Ar/rdoe.  respectively,  where  6toe  and  are  the  respective  decorre* 
ladon  parameters^  Table  1  for  N  ^(e.g.,  for  the  clean  Imaging 
system  et«e  ■©/ V2  fi  p  and  rd~  ■  2V2  /i/to). 
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TsMel.  The  Magnitude  of  the  ConeiatiODCoefiiaegl  as  a  FoBCtioa  of  die  Nundiercrf' Modes,  N. 
fiar  Three  Types  of  LoacHnrtnrwlDeftinnations 


In  contrast  to  tilt  in-plane  displacement,  jrj  decorrelates  much  slower  as  a 
function  of  Az  and  N.  These  features  are  illustrated  in  Fig.  4  where  |y|  is 
plotted  as  a  function  of  Az/z^  and  various  values  of  N.  In  particular,  for 
N««,  one  obtains  a  finite  residual  correlation  tail  equal  to  2N  /(N+1). 
These  features  are  illustrated  in  Fig.  4  where  |3r|  for  a  clean  imaging  system 
is  plotted  as  a  function  of  Az/z  ,  where  z  »  Aff/kc^  is  the  depth  of  field 
of  the  objective  lens  and  various  values  of  the  number  of  modes.  For  finite 
N,  decorrelation  effects  become  important  when  Az  >  z^.  As  an  example 
consider  wavelengths  in  the  visible  or  near  Ir  and  Az  x*  0.1  mm.  Hence,  if 
f a  10  load  induced  ‘  deformations  along  the  optic  ■  axis  are  manifest  in 
phase  difference  interferometry.  Such  a  situation  may  be  obtained,  for 
example,  in  a  system  that  employs  a  microscope. 

For  the  Fourier  transform  system  the  characteristic  decorreiation  length  is 
given  by  the  Rayleigh  range  of  the  Fourier  transform  lens.  For  practical 
applications  kr’  is  generally  much  greater  than  all  load  deformations  Az  of 
Interest  and  hence  we  conclude  that  Fourier  transform  systems  will  be 
Iniensttlve  to  displacements  along  the  optic  axis  in  phase  difference 
interferometry  applications. 
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Figure  4.  The  magnitude  of  the  correlation  coefficient  for  a  "clean” 
Imaging  system  as  a  function  of  Az/zq  for  various  values  of  N,  where 
Dz  is  the  deformation  along  the  z-axls  and  zq  ■  4/i  /  k<T*  Is  the  depth 
of  field. 
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Appendix 


Here  we  demonitrate  that  the  reiults  of  Ref.  5  for  the  correlation  function 
of  the  field  are  valid  for  N  »  1,  where  N  is  the  number  of  modes. 

Owner-Petersen  considers  a  direct  Imaie-eollectlni  geometry  where  a  lens  of 
diameter  D  collects  light  reflected  of  an  object  of  amplitude  reflection 
coefficient  p(r).  Central  tc  the  derivation  of  the  correlation  function  In 
Ref.  5  is  the  quantity  which  is  obtained  In  the  calculation  of  the 

wavefront  in  the  pupil  plane.  This  quantity  is  the  Fourier  transform  of  p  and 
in  the  notation  employed  here  is  given  by 


where 


Ap(£)  •  Jd*r  p(r)  e  ‘ 


(A-l) 

(A-2) 


L  is  the  distance  between  the  object  and  the  lens,  and  r^  is  a 

-p 

two-dimensional  coordinate  vector  in  the  pupil  plane  In  computing  the 

correlation  function  Owner-Petersen  asserts,  in  Eq.  AlO  of  Ref.  5  that  A.  is 

P 

delta  correlated: 

<Ap(aj)  6(3^-  3^)  {A-3) 

To  see  under  what  condition  this  is  valid  we  have  from  Eq.  A-l  that 


^Aptij)  Ap(aj)>  -  Jd'rJdVj  <plrj)  p»(r^)>  (A-4) 

For  a  completely  Incoherent  source  the  phase  of  p  is  random  with  statistics 
<p(r  )p*(r  )>  ■  constant  d(r  -r  J.  Using  this  into  Eq.  (A-4)  yields 

^Ap^Sj)  A*  (3j)>  ■  constant  Jd*r|p(r)|*  e*^ai“  (A-5) 
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Now  In  order  that  the  right  hand  aide  of  Eq.  A-5  resemble  a  delta  function  it 
is  necessary  for  the  exponential  term  in  the  integral  to  undergo  many 
oscillations  for  q^«  q^.  This  condition  will  bo  satisfied  if  |p|  is  siowly 
varying  and 


kDr 


max 


»  1. 


(A-6) 


where  r^  is  of  the  order  the  size  of  the  object.  Comparing  Eg.  A-6  with 
Eq.  47  indicates  that  Eq.  A-6  is  equivalent  to  the  condition  that  the  number 
of  modes  N  »  1, 


Indeed,  Eq.  A-6  can  be  obtained  by  an  elementary  argument.  Physically,  one 
can  consider  the  optical  field  in  the  pupil  plane  to  be  delta  correlated  when 
the  lateral  coherence  length  of  the  field  is  much  less  than  the  pupil 
diameter  D.  By  the  van-Clttert  Zornlke  theorem,  the  coherence  length  of  an 
incoherent  radiator  of  size  r^  is  of  the  order  Ar^/L,  where  A  it  the  optical 
wavelength.  The  condition  that  this  coherence  length  is  much  less  than  D  is 
identical  to  Eq.  A-6.  For  many  application  Eq.  A-6  la  not  satisfied  and  the 
assumption  that  A^  is  delta  correlated  is  not  valid. 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Coiporation  functions  as  an  "architect-engineer"  for  national  security 
programs.  speciaUalng  in  advanced  military  space  systems.  The  Corporation's  Technology 
Opmtlons  supports  the  effective  and  timely  development  and  operation  of  naUonal  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide*ranglng  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  Issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  Technology 
Centers: 

BUctronics  Technology  Centert  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening.  tUta  storage 
technologies,  infrwed  detector  devices  and  testing;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optical  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  solar  cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Centert  Evaluation  and  characterization  of 
new  materials:  metals,  alloys,  ceramics,  polymers  and  their  composites,  and  new 
forms  of  carbon;  development  and  analysis  of  thin  films  and  deposition  techniques; 
nondestructive  evaluation,  componmt  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  develqpment  and  evaluation  of  hardened  components; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch 
vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight  dynamics;  chemical  and 
electric  propulsion;  spacecraft  structural  mechanics,  spacecraft  survivability  and 
vulnerability  assessment;  contamination,  thermal  and  structural  control;  high 
temperature  thermomechanics,  gas  kinetics  and  radiation;  lubrication  and  surface 
phenomena. 

Space  and  Environment  Technology  Centert  Magnetospherlc,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospherlc  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  inftared  signature  analysts;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-fleld- 
of-vlew  rejection. 


